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ABSTRACT 

We analyze the gas kinematics of damped Lya systems (DLAs) hosting high z gamma-ray bursts 
(GRBs) and those toward quasars (QSO-DLAs) focusing on three statistics: (1) Awgo, the velocity 
interval encompassing 90% of the total optical depth, (2,3) IV1526 and Wcw, the rest equivalent widths 
of the Si II 1526 and C IV 1548 transitions. The Av 90 distributions of the GRB-DLAs and QSO-DLAs 
are similar, each has median Avqo rs 80 km s~ 4 and a significant tail to several hundred km s~ x . This 
suggests comparable galaxy masses for the parent populations of GRB-DLAs and QSO-DLAs and we 
infer the average dark matter halo mass of GRB galaxies is < 10 12 Mq. The unique configuration of 
GRB-DLA sightlines and the presence (and absence) of fine-structure absorption together give special 
insight into the nature of high z, protogalactic velocity fields. The data support a scenario where the 
Augo statistic reflects dynamics in the interstellar medium (ISM) and IV1526 traces motions outside 
the ISM (e.g. halo gas, galactic-scale winds). The IF1526 statistic and gas metallicity [M/H] are tightly 
correlated, especially for the QSO-DLAs: [M/H] = a + &log(W 1526 /lA) with a = -0.92 ± 0.05 and 
b = 1.41 ± 0.10. We argue that the IV1526 statistic primarily tracks dynamical motions in the halos of 
high z galaxies and interpret this correlation as a mass-metallicity relation with very similar slope to 
the trend observed in local, low-metallicity galaxies. Finally, the GRB-DLAs exhibit systematically 
larger IV1526 values (> 0.5A) than the QSO-DLAs (< IV1526 >~ 0.5A) which may suggest galactic- 
scale outflows contribute to the largest observed velocity fields. 

Subject headings: quasars : absorption lines 



1. INTRODUCTION 

Absorption-line spectra of quasars have revealed thou- 
sands of high redshift galaxies to date (e.g. Wolfe et 
al. 1986; Steidel k Sargent 1992; Prochaska, Herbert- 
Fort & Wolfe 2005; Prochter, Prochaska & Buries 2006). 
At z > 2, these galaxies are termed the damped Lya 
systems (DLAs), absorbers with H I column density 
Vhi > 2 x 10 20 cm~ 2 . Although a direct association 
between star forming galaxies and DLAs is difficult to 
establish because the background quasar blinds our view 
(Moller et al. 2002), the large N m values of DLAs in- 
dicate large dark matter overdensities (i.e. a virialized 
system). Furthermore, all of the DLAs show substan- 
tial enrichment from heavy metals (Pettini et al. 1994; 
Prochaska et al. 2003) and roughly half show indica- 
tions of ongoing star formation through the presence of 
C II* absorption (Wolfe, Prochaska & Gawiser 2003). 
Finally, stars form from neutral gas and the DLAs domi- 
nate the atomic hydrogen reservoir at all redshifts (Wolfe 
et al. 1995; Prochaska, Herbert-Fort & Wolfe 2005; Rao, 
Turnshek & Nestor 2006). 

High resolution spectroscopy of DLAs yield precise 
measurements of the column densities of resonance-line 
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transitions in the ISM of high z galaxies (Wolfe, Ga- 
wiser & Prochaska 2005). Aside from gas-phase abun- 
dances, however, it is difficult to derive physical quanti- 
ties from these observations (e.g. temperature, density). 
The only other physical characteristic easily studied is 
the gas kinematics. High-resolution spectra resolve the 
line profiles of these transitions into 'clouds' which trace 
the velocity fields in young, high z galaxies (Prochaska 
& Wolfe 1997, 1998). Although it is difficult to reveal 
the nature of these velocity fields (e.g. rotation, outflows, 
turbulence) with individual 1-dimcnsional sightlines, the 
distributions of kinematic characteristics provide power- 
ful tests for scenarios of galaxy formation (e.g. Jedamzik 
& Prochaska 1998; Haehnelt, Steinmetz & Rauch 1998). 

Prochaska & Wolfe (1997) performed a survey of the 
kinematic characteristics of neutral gas in DLAs and 
compared their observations against a variety of simple 
scenarios. These included rotating disks, clouds with 
random (virialized) motions, and galactic infall. Ob- 
served asymmetries in the observed line-profiles favor ro- 
tational dynamics, yet the observed distribution of veloc- 
ity widths (Aw) has a median too large to be accommo- 
dated within standard cold dark matter cosmology. As 
such, Prochaska & Wolfe (1997) presented these observa- 
tions as a direct challenge to the basic picture of galaxy 
formation in a hierarchical universe (e.g. Kauffmann 
1996; Mo, Mao & White 1998). In response, theorists 
introduced additional velocity fields to explain the obser- 
vations (Haehnelt, Steinmetz & Rauch 1998; McDonald 
& Miralda-Escude 1999; Nulsen, Barcons & Fabian 1998; 
Mailer et al. 2001). Of particular interest was a model of 
merging protogalactic clumps where the observed kine- 
matics include contributions from rotation, infall, and 
random motions (Haehnelt, Steinmetz & Rauch 1998). 
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While the initial work established this model as a poten- 
tial solution, recent evaluations using modern numerical 
simulations and more accurate treatments of radiative 
transfer have not confirmed its viability (Prochaska & 
Wolfe 2001; Razoumov et al. 2006). The research on 
kinematic characteristics in DLAs toward quasars (here- 
after QSO-DLAs), therefore, tests theories of galaxy for- 
mation and gas dynamics in the early universe. 

Motivated by these results, we have investigated the 
dynamics of gas probed by the sightlines to the after- 
glows of gamma-ray bursts (GRBs). Similar to quasars, 
the afterglows of GRBs provide bright, albeit transient 
point-source beacons to the outer universe. With follow- 
up spectroscopy, one can acquire data at comparable 
signal-to-noise (S/N) and spectral resolution as quasars 
(Vreeswijk et al. 2004; Fiore et al. 2005; Chen et al. 
2005; Prochaska et al. 2007). In turn, one can study 
the velocity fields of the gas close to the observed bea- 
con. The sightline presumably intersects gas in the star 
forming region encompassing the GRB, the neutral ISM 
surrounding it, and gas external to the neutral ISM. 

In this paper, we examine the kinematics of a modest 
sample of damped Lya systems associated with the ISM 
surrounding gamma-ray bursts (GRB-DLAs). A princi- 
pal goal is to characterize the velocity fields of gas ob- 
served along these unique sightlines. The data trace ve- 
locity fields generated by galactic rotation, gravitational 
accretion, and any galactic or stellar feedback processes. 
We will find that these observations reveal new insight 
for interpretations of the kinematics of gas observed in 
quasar absorption line studies. Furthermore, we will 
present a comparison of the GRB-DLA results with those 
from QSO-DLAs. Comparisons between GRB-DLAs and 
QSO-DLAs assess the relative galaxy masses and the na- 
ture of velocity fields in young, star-forming galaxies. 

This paper is organized as follows. Section 2 presents 
the observational samples and describes the experiment. 
In § 3, we define statistics for assessing the gas kinemat- 
ics. § 4 presents the principal results of our analysis and 
we discuss and interpret these results in § 5. 

2. OBSERVATIONAL DATA 

We have defined a sample of GRB-DLAs for kinematic 
analysis based on two primary criteria: (i) the pres- 
ence of a damped Lya system (Ahi > 2 x 10 20 cm~ 2 ) 
or a low-ion column density that requires TVhi exceed 
2 x 10 20 cm~ 2 assuming solar metallicity; (ii) spectra 
with sufficient signal-to-noise and resolution to study 
the gas kinematics. The latter criterion is not espe- 
cially strict because we will examine the gas kinematics 
from equivalent width measurements which do not re- 
quire high spectral resolution. The sample is also limited 
to GRB-DLAs where we could access the data or where 
precise equivalent width measurements are reported in 
the literature. Figure 1 shows an example set of line- 
profiles for the GRB-DLA associated with GRB 050820 
(Prochaska et al. 2007). These data were acquired 
with the HIRES spectrometer on the Keck telescope and 
have high spectral resolution (R sa 40,000). The weaker 
transitions such as Zn II 2026 indicate that the major- 
ity (90%) of the gas is localized to one or two 'clouds' 
within a velocity interval of ~ 50 km s _1 . In contrast, 
the stronger transitions are typically saturated and show 
absorption spanning nearly 400 km s _1 . In the following 
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Fig. 1. — Resonance-line profiles for the GRB-DLA associated 
with GRB 050820 (Prochaska et al. 2007). The velocity v = 
corresponds to z = 2.61469. Note that the optical depth in low-ion 
gas is dominated by a few clouds with velocity interval Atigo ~ 
50 km s" 1 (Zn II 2026). At the same time, there is absorption 
traced by strong low-ion transitions (e.g. Si II 1526) to velocities 
of » ±200 kms- 1 . 

section, we introduce two kinematic statistics to charac- 
terize these properties. Table 1 lists the GRB-DLA com- 
prising our sample, describes the spectral observations, 
and lists the Nm values and metallicities as described in 
Prochaska, Chen & Bloom (2007). 

Our analysis will frequently draw comparisons between 
GRB-DLA measurements and QSO-DLA values. For the 
latter, we will consider two samples of QSO-DLAs, one 
a subset of the other. The complete sample is drawn 
from echelle and echellette observations of quasars ac- 
quired with the HIRES (Vogt et al. 1994) and ESI 
(Sheinis et al. 2000) spectrometers at the Keck Obser- 
vatory and the UVES (Dekker et al. 2000) spectrometer 
at the VLT Observatory. The data are summarized in 
these papers: Herbert-Fort et al. (2006); Ledoux et al. 
(2006); Dessauges-Zavadsky et al. (2006); Prochaska et 
al. (2007). We will restrict the samples to zdla > 1.6 
and DLAs which are greater than 3000 km s _1 from their 
background quasar. The full dataset is a heterogeneous 
sample of QSO-DLAs including systems selected on the 
basis of strong metal-lines (Herbert-Fort et al. 2006) 
or targeted for H2 absorption (Ledoux, Petitjean & Sri- 
anand 2003). As such, Prochaska et al. (2007) have 
noted that the H I distribution f(Nm) of this full sample 
does not follow the statistical distribution 6 derived from 
a random survey of background quasars (e.g. Prochaska, 

6 The QSO-DLAs observed at high spectral resolution have sys- 
tematically fewer systems with TVhi ~ 2 x lO 20 cm -2 , presumably 
because the observers wished to avoid including absorbers right at 
the threshold. 
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TABLE 1 
GRB-DLA SAMPLE 



GRB 


RA 
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3.372 


21 90+ - 07 


12 
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VT T /FC)liS9 

V 1j1 / V KJIXO^i 


1 000 


7 


GRB050401 


16:31:28.82 


+02:11:14.8 


2.899 


99 fin +0 ' 30 

ZZ.DU_q 3() 


12 


-1.57 




VLT/FORS2 


1,000 


8 


GRB050505 


09:27:03.20 


+30:16:21.5 


4.275 


22 05+ - 10 


11 


-1.25 




Keck/LRIS 


1,000 


9 


GRB050730 


14:08:17.14 


-03:46:17.8 


3.969 


22 15+ - 10 
zz - lo -0 .10 


4 


-2.26 


0.14 


Magellan/MIKE 


30,000 


10 


GRB050820 


22:29:38.11 


+19:33:37.1 


2.615 


21 00+ - 10 
zi.uu_ .i 


4 


-0.63 


0.11 


Keck/HIRES 


30,000 


11 


GRB050904 


00:54:50.79 


+14:05:09.4 


6.296 


21 30+ - 20 
il - ou -0 .20 


11 


-1.10 




Subaru/FOCAS 


1,000 


12 


GRB050922C 


19:55:54.48 


-08:45:27.5 


2.199 


21 60+ - 10 


4 


-2.03 


0.14 


VLT/UVES 


30,000 


13 


GRB051111 


00:08:17.14 


-00:46:17.8 


1.549 


22 6 


12 


-0.96 




Keck/HIRES 


30,000 


14,11 


GRB060206 


13:31:43.42 


+35:03:03.6 


4.048 




4 


-0.85 


0.18 


WHT/ISIS 


4,000 


15 


GRB060418 


15:45:42.40 


-03:38:22.80 


1.490 


22 6 


2 


-1.65 


1.00 


Magellan/MIKE 


30,000 


11 



REFERENCES. — 1: Savaglio, Fall & Fiore (2003); 2: Castro ei al. (2003); 3: Mirabal ei al. (2002); 4: Vrceswijk et al. (2006); 5: Barth et al. (2003); 6: Shin 
et oi. (2006); 7: Vrccswijk et a/. (2004); 8: Watson et a/. (2006); 9: Bcrgcr et a/. (2006); 10: Chen et al. (2005); 11: Prochaska et al. (2007); 12: Kawai at al. 
(2006); 13: Piranomontc et al. (2007); 14: Prochaska, Chen & Bloom (2006); 15: Fynbo at al. (2006) 

a Flag describing the mctallicity measurement [M/H]: = No measurement; l = Si measurement; 2 = Zn measurement; 3 = Combination of limits; 4=S measurement; 
ll=Lower limit from [a/H]; 12 = Lowcr limit from Zn 

b Bceausc z GRB < 1.6, Lyct was not observed. We have set jV m = 10 22 cm -2 , the median value of GRB-DLA 



Herbert-Fort & Wolfe 2005). Therefore, we define a 
pseudo-statistical sample of QSO-DLAs by restricting 
the list to the QSO-DLAs compiled by Prochaska et al. 
(2003). Although the DLA in this subset also do not 
follow the /(./Vhi) distribution of a random sample, they 
were selected only on the basis of a large H I column 
density, i.e. independent of any kinematic characteristic 
or chemical abundance. 

We have analyzed the kinematic properties from our 
own HIRES and ESI observations and supplemented 
these measurements with the results presented in Ledoux 
et al. (2006) from UVES data acquired with the VLT. 
For equivalent width measurements, however, we have 
derived values only from our HIRES and ESI observa- 
tions (Prochaska et al. 2007; Herbert- Fort et al. 2006). 
These values have typical statistical uncertainties of less 
than 0.02A. 

In a previous paper (Prochaska, Chen & Bloom 2007) , 
we discussed why the parent populations of galaxies 
hosting QSO-DLAs and GRB-DLAs may be different 
(e.g. gas cross-section selected versus current SFR se- 
lected) and how the sightline configurations are funda- 
mentally distinct. The GRB-DLAs do exhibit larger -ZVhi, 
metallicity, a/Fe, and depletion levels than QSO-DLAs 
(Prochaska, Chen & Bloom 2007), but we argued that 
the differences do not require distinct parent populations 
of galaxies. Instead, the observations imply sightlines 
with smaller average impact parameter through high z 
galaxies. This conclusion is expected for a sample se- 
lected by gas cross-section (QSO-DLA) versus ones re- 
stricted to originate in star- forming regions (GRB-DLA). 
The cartoon in Figure 2 illustrates these ideas. While 
QSO-DLAs preferentially penetrate the outer regions of 
the ISM, the GRB-DLA sightlines originate within their 
host galaxies, presumably within an H II region gener- 
ated by its progenitor and other O and B stars. This 




Fig. 2. — This cartoon illustrates the likely differences between 
QSO-DLA and GRB-DLA sightlines. The former have randomly 
intersected a foreground galaxy. These QSO-DLA sightlines corre- 
spond to a cross-section selected sample and should preferentially 
intersect the outer regions of the ISM in high z galaxies. In con- 
trast, the GRB-DLA sightlines are constrained to originate from 
within the ISM of their host galaxies, presumably the H II region 
produced by massive stars in a star-forming region. These GRB- 
DLA sightlines are expected (and observed) to originate within the 
inner few kpc of the ISM. 

has two key implications for the gas kinematics: (i) the 
GRB-DLA observations will probe only a fraction of the 
velocity field along the full sightline; (ii) each GRB-DLA 
sightline is biased to intersect its own star forming region. 
We expect these regions to have very small cross-section 
to QSO-DLA sightlines (Zwaan & Prochaska 2006) and 
to be rarely probed by such samples. In this respect, 
GRB-DLAs may open a new window into the velocity 
fields of high z galaxies (c.f. Pettini et al. 2001). 

3. KINEMATIC DIAGNOSTICS 

With high resolution observations, one resolves the 
absorption-line profiles of unsaturated transitions into 
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TABLE 2 

GRB-DLA KINEMATIC CHARACTERISTICS SUMMARY 



GRB000926 


2026.136 


280 


2.40 ±0.14 


> 2.20 


GRB010222 


2056.254 


100 


1.51 ±0.13 


> 2.55 


GRB011211 






1.35 ±0.19 


0.70 ± 0.13 


GRB020813 


2260.780 


210 


1.71 ±0.10 


1.53 ±0.05 


GRB030226 


2374.461 


340 a 


0.97 ± 0.02 


0.32 ± 0.02 


GRB030323 






0.75 ± 0.03 


1.50 ±0.03 


GRB050401 






2.31 ±0.26 


1.28 ±0.26 


GRB050505 






1.81 ±0.10 


> 5.92 


GRB050730 


1741.553 


25 


0.37 ±0.01 


0.81 ±0.01 


GRB050820 


2026.136 


55 


1.65 ±0.01 


1.51 ±0.01 


GRB050922C 


1608.451 


89 


0.52 ±0.01 


0.75 ±0.01 


GRB051111 


2249.877 


31 


0.79 ± 0.20 


1.35 ±0.20 


GRB060418 


2576.877 


57 


0.66 ± 0.02 


0.80 ±0.02 


Note. — All EW v 






equivalent width 




a Estimatcd from F 


gurc 2 of Shii 


i et al. (2006). 


Because the lin 




tially saturated, this 


value may be 


considered an 


upper limit. 





individual components, frequently termed 'clouds'. The 
clouds in QSO-DLAs have typical Doppler widths of 
b = \/2<t sa 5 — 10 km s _1 (e.g. Dessauges-Zavadsky et al. 
2006). The dispersion results from macroscopic motions 
(e.g. turbulence) not thermal broadening; the latter im- 
ply temperatures that would collisionally ionize the gas. 
The GRB-DLAs show clouds with similar characteristics 
that comprise the observed line- profiles (Figure 1; Fiore 
et al. 2005; Chen et al. 2005; Prochaska et al. 2007). 

One can characterize the velocity field along the sight- 
line by tracing the motions of these clouds. Presently, we 
do not have measurements of the systemic redshifts for 
our sample of GRB host galaxies (e.g. via nebular lines). 
Therefore, we will primarily discuss relative velocities. 
Standard practice is to define a velocity width Avqq as 
the interval which encompasses 90% of the total optical 
depth of the gas. This definition was introduced in part 
to limit a single, weak cloud or statistical fluctuations 
from dominating the statistic (Prochaska & Wolfe 1997). 
Following the approach for QSO-DLAs (Prochaska & 
Wolfe 1997, 1998, 2001; Ledoux et al. 2006), we mea- 
sure the Avgo statistic from a single, unsaturated low- 
ion transition 7 which, ideally, has relatively high S/N. 
For the QSO-DLA we demand S/N > ISpix" 1 but have 
relaxed this requirement for the GRB-DLA because of 
the generally poorer data quality. To resolve the optical- 
depth profile, however, it is necessary to restrict the sam- 
ple to echelle or echellette observations. As an example of 
the procedure, consider the line-profiles in Figure 1 where 
we measure Avgo = 55 km s^ 1 using the Zn II A2026 
transition. Table 2 summarizes the measurements for 
the full GRB-DLA sample. 

We also introduce a complimentary diagnostic of gas 
kinematics: the rest equivalent width W = W b s / (1 + z). 
This quantity has kinematic significance for optically 
thick (i.e. saturated) transitions. In these cases, the 
equivalent width is most sensitive to the differential mo- 
tions of individual clouds as opposed to their column 
densities. In contrast to Afgrj, the W value represents 

7 Note that Augo is insensitive to the specific transition analyzed 
provided the line is unsaturated and is associated with a low-ion 
(Prochaska & Wolfe 1997). 
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Fig. 3. — Curve-of-growth for the Si II 1526 transition assuming 
a Doppler parameter 6 = 7 km a -1 . Overplotted on the figure are 
the observed total Si+ column densities and equivalent widths from 
the QSO-DLAs (green, plus-signs) and GRB-DLAs (blue stars). 
Aside from the few QSO-DLAs with very low W1526 , the equivalent 
widths arc dominated by the contribution from multiple clouds 
with a broad range of column densities. 

the velocity field of greater than 90% of the gas and can 
be dominated by clouds with relatively low column den- 
sity (e.g. Bouche et al. 2007). We emphasize that the W 
value should only be sensitive to the metallicity and/or 
Ahi value of the sightline at very low values, i.e. when 
the transition is not saturated. Below, we demonstrate 
that Avgo and W are correlated but with large scatter. 

We have measured rest equivalent widths from the 
Si II A1526 transition (Wi^o) for low-ion gas and the 
C IV 1548 transition for high-ion gas. The Si II A1526 
transition saturates at ./V(Si + ) w 10 14 cm~ 2 which is 
10 x lower than the Si + column density for any GRB- 
DLA sightline and 2x lower than nearly every QSO- 
DLA. Therefore, this measure should be relatively in- 
sensitive to the gas metallicity or A^ji value. This point 
is emphasized in Figure 3 where we show the curve-of- 
growth for the Si II 1526 transition and the observed 
total Si + column density and equivalent widths for the 
QSO-DLAs and GRB-DLAs. It is evident that the equiv- 
alent width has significant contributions from multiple, 
lower column density clouds. Because C IV is a doublet, 
the C IV 1548 transition will blend with C IV 1550 when 
^1548 ^ 2 A. In these few cases, we adopt a lower limit 
for Wciv- The W 152 q and W C w values of the GRB-DLA 
are summarized in Table 2. We also list the metallicity 
[M/H] measurement (or lower limit) for each GRB-DLA 
(see Prochaska, Chen & Bloom 2007, for details). 

4. RESULTS 

In this section, we present the observational results. 
The next section will discuss the implications. 

4.1. Awgo and VF1526 Distributions 

Figure 4a presents the velocity width (Ai>go) distribu- 
tions for the GRB-DLAs and for the statistical subset 
of QSO-DLAs. The two distributions are similar; the 
majority of galaxies have Avqq < 100 km s _1 and each 
sample shows a significant tail out to several hundred 
km s . A two-sided Kolmogorov-Smirnov test rules out 
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Fig. 4. — Histograms of low-ion kinematic statistics for the GRB- 
DLAs (dark blue) compared against a statistical sample of QSO- 
DLAs (light green). The upper panel shows the results for the Avgo 
statistic, the velocity interval containing 90% of the optical depth 
of the low-ion gas. The two samples have similar distributions 
with median A«go ~ 80 km s — 1 and a significant tail to several 
hundred km s — x . The lower panel shows the equivalent widths 
of the Si II 1526 transition. In contrast to the Avgo statistic, 
the GRB-DLAs exhibit systematically larger W1526 values than 
the QSO-DLAs. In both panels, the inset shows the cumulative 
distributions of the systems, normalized to unity, and the Pks 
value gives the probability that the two distributions are drawn 
from the same parent population. 

the null hypothesis that the two distributions are drawn 
from the same parent population at only 32% c.l. In 
the current GRB-DLA sample, there is a higher fraction 
of systems with Avgn > 200 km s -1 than the QSO-DLA 
distribution. The difference is not statistically significant 
at present, but future observations may reveal a system- 
atic difference in this respect. 

Because the GRB are embedded within the interstel- 
lar medium (Figure 2), the velocity widths are likely an 
underestimate of the value that one would derive by ex- 
tending the sightline through the other side of the galaxy. 
If the kinematics are dominated by random motions or 
rotational dynamics, then the median increase in Augn 
should be less than a factor of two. To be conserva- 
tive, we have recomputed the K-S statistic after dou- 
bling each Avgo value of the GRB-DLAs. We find re- 
sults that are still consistent with the null hypothesis, 
P KS = 0.08. In this regard, the sample of GRB-DLA 
sightlines have velocity fields typical of those observed 
for QSO-DLAs. Under the assumption that Avgo is 
dominated by the gravitational potential (e.g. Haehnclt, 
Steinmetz & Rauch 1998), these results suggest that the 
parent populations of GRB-DLA and QSO-DLA have 
similar dynamical masses. 

From even the first afterglow spectrum (Metzger et al. 



1997), it was evident that the gas surrounding GRB ex- 
hibit large equivalent widths from low-ion transitions like 
Si II A1526 and Mg II A2796 (Savaglio, Fall & Fiore 
2003). As noted in § 3, these large equivalent widths 
do not require large Ajji or metallicity; the values are 
often dominated by the velocity fields of low column den- 
sity clouds which do not contribute to the Avgo statis- 
tic. Figure 4b presents a histogram of \V152e values for 
the GRB-DLA compared against the values of the sta- 
tistical QSO-DLA subset. The mean and median of the 
GRB-DLA distributions are significantly larger than the 
QSO-DLAs and a K-S test rules out the null hypothe- 
sis at 99.9%c.l. Whereas the Avgn values are in rather 
good agreement between the two samples, the W1526 val- 
ues of the GRB-DLAs are systematically larger. This is 
an unexpected and puzzling result which requires unique 
velocity fields from low column density gas along GRB 
sightlines. 
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Fig. 5.— Plot of Av 90 ,W 152 6 pairs for the QSO-DLAs (light 
green) and GRB-DLAs (dark blue) samples. The red dashed-linc 
indicates the W1526 value for a boxcar line-profile with velocity 
width AD90, i.e. W1526 = Augo * 1526.7066/c. The sightlines that 
lie significantly above this curve must have contributions to W1526 
from gas at large velocity and with less than 10% of the optical 
depth. 

We will argue below that the Ai; go statistic in GRB- 
DLAs is dominated by the ISM of the GRB host galaxy. 
It is our expectation that this also holds true for QSO- 
DLAs. We will also argue that the W1526 values result 
from motions independent of the neutral ISM. If other en- 
vironments contribute to W1526, then one predicts Ai;go 
to be only loosely correlated with W^26- Figure 5, which 
plots Aiig ,Wi526 pairs for the GRB-DLAs and QSO- 
DLAs, indicates that this is the case. The dashed line 
traces the predicted W1526 value for a fully saturated 
(i.e. boxcar) line-profile that has velocity width Avqq. 
Although the QSO-DLAs roughly follow this line, the 
GRB-DLAs tend to lie significantly above it. Regarding 
the QSO-DLAs, we interpret the results in Figure 5 as 
evidence that many sightlines penetrate gas with velocity 
fields that are distinct from the majority of gas. Further- 
more, this fraction appears to increase with Aw go . In the 
GRB-DLAs, nearly every sightline shows a substantial 
contribution to W1526 from gas that does not dominate 
the total optical depth. In § 5, we will discuss the origin 
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Fig. 6. — Plots indicating correlations between the gas kinematics and metallicity. Upper panel shows pairs of Avgo, [M/H] values 
for the QSO-DLAs (light green plus-signs) and GRB-DLAs (dark blue stars). The GRB-DLA with lower limits to [M/H] have saturated 
metal-line profiles while the points marked with squares do not have Nm measurements because they have z a i, a < 1.6. The QSO-DLAs 
exhibit a previously discussed correlation (Wolfe & Prochaska 1998; Ledoux et al. 2006) and the GRB-DLA appear to trace the same 
locus. The lower panel compares [M/H] with W\s26- The QSO-DLA exhibit a remarkably tight trend, described by the dashed red line: 
[M/H] = —0.92 + 1.41 log(Fyi526/lA). The GRB-DLAs also appear to show a correlation which may be offset or steeper than the QSO-DLA, 
but this conclusion is complicated by the preponderance of lower limits. Finally, the black diamond shows the [M/H],H / i526 values for 
the Lyman break galaxy CB-58 (Pettini et al. 2002). Although its correspondence with the QSO-DLA trend may be coincidence, it is 
suggestive that large W1526 values are related to galactic outflows. 



and nature of these velocity fields in greater detail. 

4.2. Kinematic Correlations with Other Properties of 

the ISM 

We have considered the trends between kinematic char- 
acteristics and other physical properties of the ISM. We 
first examined whether the GRB-DLAs follow the ob- 
served trends for QSO-DLAs between gas kinematics and 
metallicity (Wolfe & Prochaska 1998; Ledoux et al. 2006; 
Murphy et al. 2007). Figure 6a presents [M/H] values 
against the Avgo statistic for the GRB-DLAs and the 
full QSO-DLA sample. Note that the GRB-DLA data 
points with lower limits to [M/H] have z a bs < 1.6 and an 
assumed H I column density A/hi — 10 22 cm~ 2 . There- 
fore, these [M/H] values (reported as lower limits because 
their metal-line profiles are saturated) should be consid- 
ered cautiously. The QSO-DLAs shown in Figure 6a ex- 
hibit the correlation reported by previous authors (Wolfe 
& Prochaska 1998; Ledoux et al. 2006). It is not a tight 
trend, however, and we suspect several physical factors 
(e.g. a mass/metallicity relation; variations with sightline 
impact-parameter and galaxy inclination) contribute to 
produce the observed distribution. Similar to the Avgo 
distribution (Figure 6a), we find that the GRB-DLAs 
track the locus defined by the QSO-DLAs. Adopting the 
upper limits to [M/H] as values, we report a Pearson 
correlation coefficient of 0.4. We conclude that there is 
only tentative evidence for a correlation between Avgo 
and [M/H] for the GRB-DLA. 

In Figure 6b, we present the W1526) [M/H] pairs for 
the two populations. The QSO-DLA data exhibit a re- 
markably tight correlation. These are are well described 
by a power-law, 

[M/H] =o + 61og(W/lA) , (1) 

with best-fit parameters, a = —0.92 ± 0.05 and b = 
1.41 ± 0.10. For this fit, we have assumed equal weights 
for all of the data points and we have restricted the 
sample to W1526 < 1.4A to ignore the 'outliers' at 
large IF1526 value. If we include all of the data points 
(with equal weighting), we derive a = —1.00 ± 0.05 and 
b = 1.27 ± 0.10. Although the data scatter about this 
power-law, the trend is considerably tighter than the 
correlation between Align and [M/H]. In fact, it is the 
tightest correlation known between metallicity and any 
other property of the QSO-DLAs (see also Murphy et al. 
2007). The small scatter is especially impressive given 
that « 0.15 dex is expected from observational uncer- 
tainty in the [M/H] measurements (A/hi error). It is a 
rather surprising trend because the kinematics of the gas 
that determines [M/H] are better described by the Avgo 
statistic yet Figure 6a shows this correlation has much 
greater scatter. Put another way, the line-profiles with 
W1526 > 0.3A are highly saturated and should not be 
expected to reflect the gas metallicity (Figure 3). 

This is not to imply that the W1526 values are strictly 
independent of [M/H]. For example, an optical depth 
profile characterized by shallow, decreasing wings ex- 
tending to large velocity (e.g. a Lorentzian profile) would 
yield larger W1526 for larger A(Si + ). On the other hand, 
a velocity field with a sharp cutoff (e.g. a Gaussian ran- 
dom field) would be insensitive to [M/H]. Furthermore, 
the peak optical depth value is a function of both [M/H] 
and A/hi, and the QSO-DLA data presented in Figure 6b 
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Fig. 7.— Pairs of W1526, N m values for the QSO-DLAs (light 
green) and GRB-DLAs (dark blue). Regarding the QSO-DLA, sys- 
tems with larger W1526 show larger Nm value on average, yet the 
correlation shows large scatter. Similarly, the GRB-DLA sample 
is characterized by a large scatter but without any obvious trend. 
Although the Afpji value and metallicity [M/H] contribute equally 
to the Si + column density and presumably W1526, the latter is 
tightly correlated only with metallicity. 



span a factor of « 1.5 dex in Nm value. Therefore, the 
results on the QSO-DLAs in Figure 6b indicate an under- 
lying physical mechanism which causally connects [M /H] 
and I4 7 i526- The GRB-DLAs also show a trend between 
W1526 and [M/H], although its characterization is com- 
plicated by the many lower limits to [M/H]. A qualitative 
assessment of the current results suggests either an off- 
set between the trends for the two populations and/or a 
steeper power-law for the GRB-DLAs. We will return to 
these issues in § 5.3. 

Irrespective of the physical origin of the correlation, 
the results presented in Figure 6b indicate one may use 
W1526 as a proxy for metallicity, especially if applied to 
a sample of data points instead of individual systems. 
We find that a sample of at least 10 QSO-DLA systems 
gives an average logarithmic metallicity derived from the 
observed W1526 values and Equation 1 that lies within 
20% of the actual value for 95% of random trials. Em- 
boldened by this result, we apply equation 1 to the GRB- 
DLAs but offsetting [M/H] by Smb to allow for an offset 
between the GRB-DLA and QSO-DLA trends. We de- 
rive an average logarithmic metallicity for the GRB of 
< [M/H] >= -0.8 + S M B- 

We have also examined trends between VF1526 and the 
H I column density (for Avgo vs. A/hi, see Wolfe & 
Prochaska 1998; Wolfe, Gawiser & Prochaska 2005). Fig- 
ure 7 presents the Wi526,Nm pairs for the QSO-DLAs 
and GRB-DLAs. The two quantities are correlated in 
that the median A/hi value increases with Wi526- Fur- 
thermore, there are no sightlines with low VF1526 and 
large Nm value. For sightlines with W1526 > 0.3A, how- 
ever, we observe a wide range of A/hi values for any 
given W1526 value. The large observed scatter contrasts 
with the [M/H],Wi526 trends in Figure 6b, even though 
A(Si+) is the product of [M/H] and N m . We conclude 
that only [M/H] is physically tied to the gas kinematics 
expressed by the W1526 statistic. 
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Fig. 8. — Gas kinematics of high-ions (Wciv) compared against 
low- ion kinematics (W152&) for the QSO-DLA (light green) and 
GRB-DLA (dark blue) samples. The two quantities are clearly 
correlated for the QSO-DLAs, whereas the GRB-DLAs exhibit too 
small of a dynamic range to robustly exhibit any trend. 

4.3. High-ion gas: Wqtv 

It is difficult to measure the 90% velocity interval for 
C IV in GRB-DLAs because the C IV doublet is often sat- 
urated. Therefore, we have examined the rest equivalent 
width of the C IV 1548 transition, Wcrv- These are plot- 
ted against W1526 in Figure 8. Even though unsaturated 
low-ion profiles do not closely track the C IV line-profiles 
in QSO-DLAs (i.e. component by component; Wolfe & 
Prochaska 2000), there is a strong correlation between 
their low and high- ion kinematics. Wolfe & Prochaska 
(2000) have interpreted this global correlation as an in- 
dication that the gas traces the same gravitational po- 
tential (see also Mailer et al. 2003). 

The GRB-DLAs also exhibit a significant correlation 
between W1526 and Wciv- But, one also finds that the 
Wciv values of the GRB-DLA are systematically larger 
than those for the QSO-DLAs. Similarly, Wciv corre- 
lates with the ISM metallicity for both samples. The 
correlation between Wciv and W1526 in the GRB-DLAs 
suggests that their velocity fields are at least coupled. 
This assertion is further supported by the coincidence 
of C IV and Si II A1526 absorption at some velocities 
(Prochaska et al. 2007). Therefore, investigations on the 
velocity field of one ion may give insight to the other. 

5. DISCUSSION 

5.1. Summary of Results 

The results presented in the previous section can be 
summarized as follows: 

1. The distributions of Awgo values (the velocity inter- 
val encompassing 90% of the gas) from QSO-DLAs 
and GRB-DLAs are very similar. Both samples ex- 
hibit median Avqq ps 80 km s _1 with a significant 
tail to several hundred km s~ x (Figure 4a). 

2. In contrast to the Augo statistic, the GRB-DLA 
exhibit systematically larger W1526 measurements 
than the QSO-DLAs (Figure 4b). 

3. Both the GRB-DLAs and QSO-DLAs show a cor- 
relation between Avqq and metallicity [M/H] but 



with large scatter (Figure 6a). 

4. There is a remarkably tight correlation between 
W1526 and metallicity [M/H] for the QSO-DLA 
sample. The correlation is well described by a 
power-law, [M/H] = a + 61og(Wi526/lA) with a = 
-0.92 ± 0.05 and b = 1.41 ± 0.10. There is also a 
correlation observed between metallicity and W1526 
for the GRB-DLAs. This trend may be offset 
and/or steeper than the QSO-DLA distribution 
(Figure 6b). 

5. The high-ion and low- ion kinematics traced by 
W1526 and Wciv are correlated (Figure 8). 

We will now explore the origin and nature of the velocity 
fields in QSO-DLAs and GRB-DLAs and the implica- 
tions for the galaxies which produce them. 

5.2. The Nature of Galactic Velocity Fields at High z 

In a star-forming galaxy, there are four principal areas 
that have unique velocity field characteristics and can 
contribute to the Awgo and W1526 statistics: (i) the neu- 
tral ISM, whose motions are dominated by galactic rota- 
tion and mild turbulent velocities; (ii) H II regions, with 
velocity fields driven by stellar processes (e.g. winds, su- 
pernovae); (iii) halo gas clouds, whose kinematics include 
contributions from gravitational accretion and virialized 
motions; and (iv) galactic-scale outflows driven by super- 
novae activity in the ISM (e.g. Mac Low & Ferrara 1999) 
and/or heating by galaxy-galaxy mergers (e.g. Cox et al. 
2006). In the following, we will examine our observa- 
tional results within the framework of a model where the 
Afgo statistic is dominated by gas in the ambient, neu- 
tral ISM while halo dynamics determine the the W1526 
and Wciv measurements. In addition, we will propose 
that the W1526 and Wciv statistics occasionally trace 
velocity fields representative of galactic-scale outflows. 

This two component scenario (ISM+halo gas) is 
neither an original nor provocative interpretation of 
absorption-line systems arising in galaxies (e.g. Mo & 
Miralda-Escude 1996; Kauffmann 1996; Haehnelt, Stein- 
mctz k Rauch 1998; McDonald & Miralda-Escude 1999; 
Mailer et al. 2001, 2003). Indeed, this is also the fa- 
vored model for gas in the Milky Way (e.g. Spitzer 1990; 
Savage, Sembach & Lu 1997) and Magellanic Clouds (e.g. 
Lehner & Howk 2007). And, the model has been invoked 
for a wide range of quasar absorption-line observations 
including the impact parameter distribution of Mg II sys- 
tems from their host galaxies (Bergeron & Boisse 1991; 
Lanzetta & Bowen 1990, 1992; Steidel 1993) and C IV 
absorption by galaxies (Chen, Lanzetta & Webb 2001; 
Wolfe & Prochaska 2000; Mailer et al. 2003). We will 
demonstrate that our observations support this model 
and we will derive new insights into the gas dynamics of 
high z galaxies. 

Within this disk/halo scenario, the Awgo statistic 
traces the kinematics of the ambient ISM, i.e., the optical 
depth of the gas external to the neutral ISM contributes 
less than 10% of the total. In the context of QSO-DLAs, 
this is reasonable because (i) the large iVm values that 
define a DLA may require the sightlincs intersect the am- 
bient ISM; and (ii) one expects halo gas to be predomi- 
nantly ionized (e.g. Mo & Miralda-Escude 1996) whereas 
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Fig. 9. — These panels show one weak line (top), one strong line (middle) and one fine-structure transition (bottom) for four GRB-DLAs 
observed at high spectral resolution: (a) GRB-DLA/050922C; (b) GRB-DLA/060418; (c) GRB-DLA/051111; and (d) GRB-DLA/050820. 
In all four GRB-DLAs, the fine-structure optical depth profile closely matches that of the weak transition. Because this gas must lie within 
1 kpc of the GRB afterglow, these observations indicate the Avgo statistic is determined by the velocity field of the ISM surrounding 
the GRB. In contrast, there are clouds in the wings of the strong transitions which do not show significant fine-structure absorption and 
must therefore lie at distance > 1 kpc from the GRB afterglow. Because the fine-structure absorption defines the systemic velocity of the 
ISM, the unique geometry of the GRB-DLA sightlincs allows one to investigate gas with velocities indicate of outflow (negative) or infall 
(positive). This experiment offers an unprecedented view into the nature of gas kinematics in high z galaxies. 
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the majority of gas in QSO-DLA is neutral (Vladilo et al. 
2001; Prochaska et al. 2002). We will now demonstrate 
that observations of the GRB-DLAs directly connect the 
neutral ISM with the Avg statistic. 

Figure 9 plots pairs of strong resonance- line transitions 
(one weak, one strong) against a fine-structure transition 
for four GRB-DLAs observed at high spectral resolution. 
In all four cases, the fine-structure absorption lines trace 
only the gas which dominates the total optical depth, 
i.e. the gas revealed by the weak transitions. Prochaska, 
Chen & Bloom (2006) showed that fine-structure absorp- 
tion in the GRB-DLAs is due to UV pumping by radi- 
ation from the GRB event and its afterglow. For the 
GRB-DLAs shown here, the presence of fine-structure 
absorption dictates that this gas is located within « 1 kpc 
of the burst. Furthermore, the GRB-DLAs also show 
Mg I absorption coincident with the fine-structure lines 
(Prochaska et al. 2007). This requires the gas to lie at 
> 100 pc to avoid being photoionized by the GRB af- 
terglow (Prochaska, Chen & Bloom 2006, 2007). The 
overall implication is that gas comprising 90% of the to- 
tal optical depth lies between w 100 pc and w 1 kpc from 
the GRB, i.e. within the nearby ISM 8 . In summation, 
the GRB-DLA observations give empirical confirmation 
that the Avg values reflect the kinematics of the neutral 
ISM. 

Having established a direct association between the 
At>go values and the ISM, it is reasonable to assume that 
these values trace the rotational dynamics and turbulent 
fields within high z galaxies. This is especially the case 
for systems with At>g < 200 km s _1 ; larger Avg values 
may be better interpreted as sightlines through multiple 
galaxies or even supernovae winds (Nulscn, Barcons & 
Fabian 1998). Because the turbulent velocities in neu- 
tral gas are generally small (a tur b ~ 10 km s _1 ), sight- 
lines with A^go > 20kms _1 should be dominated by 
differential rotational along the sightlinc or peculiar mo- 
tions of multiple massive clumps (Haehnelt, Steinmctz 
& Rauch 1998; Mailer et al. 2001). It is natural, there- 
fore, to compare the Awg distributions of GRB-DLAs 
and QSO-DLAs in terms of dynamical masses. Before 
drawing such conclusions using the Awgo distributions, 
however, one must examine the effects of differences in 
sightline configuration (Figure 2). 

The GRB-DLA sightlines are restricted to travel 
through approximately half of the galaxy and are ex- 
pected to have smaller impact parameter than the aver- 
age QSO-DLA sightline. The latter point is supported by 
the higher H I column densities, metallicities, and deple- 
tion levels observed for GRB-DLA sightlines (Prochaska, 
Chen & Bloom 2007). These two effects will compete 
against each other if differential rotation explains the 
Awgo statistic. On the one hand, a reduced sightline 
leads to an underestimate of Awgo because one does not 
probe all of the gas along the sightline. This effect 
should be small, however, because even a full sightline 
will only probe approximately one quadrant of a rotating 
disk (Prochaska & Wolfe 1997). The correction would be 
larger for a random velocity field, but still less than a fac- 
tor of two. It is only if the velocity fields are dominated 

8 The absence of H2 absorption in this gas (Tumlinson et al. 
2007) also indicates we are probing the ambient ISM and not the 
local star forming region (i.e. molecular cloud). 



by a symmetric infall/outflow that the effect would be 
large. Such velocity fields, however, are unlikely to de- 
scribe the neutral ISM of high z galaxies. The other ge- 
ometric effect - smaller average impact parameter - will 
imply larger Av go values (e.g. Prochaska & Wolfe 1997). 
This effect is likely to dominate and, if anything, we ex- 
pect the GRB-DLAs would be biased to higher A^go val- 
ues than the QSO-DLAs for sightlines penetrating the 
same parent population of galaxies. 

We conclude, therefore, that the observed A^go dis- 
tributions indicate the galaxies hosting QSO-DLAs and 
GRB-DLAs have similar mass distributions. There is at 
least circumstantial evidence in support of this conclu- 
sion. In particular, the median luminosity of GRB-DLAs 
host galaxies is low, approximately L,/5 at z ~ 1 (Lc 
Floc'h et al. 2003). Although very few QSO-DLA have 
been successfully imaged at z > 2 (M0ller et al. 2002), 
theoretical expectation is that the majority will also 
be sub-luminous (Haehnelt, Stcinmetz & Rauch 2000; 
Mailer et al. 2000; Nagamine et al. 2005) and the aver- 
age luminosity at z < 1 is L <~ L*/2 (Chen & Lanzetta 
2003). Similarly, neither sample has shown evidence 
for a large fraction of bright galaxies (Le Floc'h et al. 
2003; Colbert & Malkan 2002). Altogether, we conclude 
that the host galaxies of GRB-DLAs and QSO-DLAs 
have comparable masses. Adopting recent results on the 
latter population, we infer dark matter halo masses of 
< 10 12 M Q for the GRB-DLA (Bouche et al. 2005; Cooke 
et al. 2006; Nagamine et al. 2005). 

The assertion that GRB-DLA and QSO-DLA host 
galaxies have similar mass is immediately challenged 
by our observation that GRB-DLA show systematically 
larger values (Figure 4b). To address this issue, we 

must identify the origin of the velocity fields contributing 
to Wi526- There are two factors which contribute to the 
IL1526 statistic: the number of clouds intersected by the 
sightline and the relative velocities of these clouds. The 
total velocity field is a convolution of the relative velocity 
field within each environment giving rise to Si + and the 
relative velocities between the environments. We empha- 
size that even highly ionized gas can have trace quantities 
of Si + that would be revealed by the strong Si II A1526 
transition (e.g. Prochaska 1999). 

Returning to Figure 9, we can address the contribu- 
tion of various velocity fields to Wx^s- At small W1526 
values (a,b), the neutral ISM clearly dominates. In con- 
trast, the two GRB-DLAs with large IF1526 (c,d) show 
multiple clouds with relatively low column density and no 
associated fine-structure absorption. The absence of fine- 
structure absorption places these clouds at > 1 kpc from 
the GRB afterglow (Prochaska, Chen & Bloom 2006; 
Chen et al. 2007). Note that this is true even though the 
column densities of these components are small compared 
to the total because the resonance and fine-strucure tran- 
sitions have similar oscillator strengths. Therefore, this 
gas cannot be associated with star-formation processes 
local to the GRB (e.g. stellar winds, supernovae). The 
gas is either within the ISM but at much larger distance 
than the gas that dominates the total optical depth, or 
it lies outside the ISM altogether (i.e. halo gas or the 
ISM of a companion galaxy). We prefer to associate the 
gas with the galactic halo because: (i) the velocities are 
large (up to hundreds km s _1 ) compared to random and 
dynamical motions for the neutral ISM of galaxies; (ii) 
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if the W1526 statistic tracked dynamics in the ISM then 
a very large metallicity gradient may be required for the 
Ai; 90 statistic to be insensitive to these motions; (hi) 
there is high-ion absorption (e.g. Si IV, C IV; Figure 1) 
coincident with many of the low column density Si + com- 
ponents; and (iv) the ionic ratios of some components 
indicate they are highly ionized (Dessauges-Zavadsky et 
al. 2006). An example of the last point is revealed in 
Figure 1 by the cloud at v « +200 km s _1 . This cloud 
has log[iV(O )/iV(Si + )] = +0.50 ±0.05 even though O is 
at least 1 dex more abundant than Si in every astrophysi- 
cal environment known. The only viable explanation for 
observing O /Si + <C (0/Si)o is that the gas is highly 
ionized. 

Granted that gas beyond the ISM dominates the W1526 
statistic, the next issue is to examine the nature of its ve- 
locity field. This may include virializcd motions, gravita- 
tional accretion (Mo & Miralda-Escude 1996; McDonald 
& Miralda-Escude 1999), and outflows from galactic feed- 
back (Dong, Lin & Murray 2003; Cox et al. 2006). The 
latter two processes imply organized flows that would 
be revealed by the observations if we knew the systemic 
velocity of the galaxy. For the GRB-DLAs, we can mea- 
sure the systemic redshift from nebular emission lines 
(e.g. Thoene et al. 2006) but the absorption-line data 
also yield a direct measurement: the line centroid of the 
fine-structure transitions. Because this gas occurs within 
« 1 kpc for the GRB-DLA, it establishes the velocity of 
the neutral ISM. Examining the line-profiles of Figure 9 
(c,d) in this light, we note that the Si + gas shows both 
positive and negative absorption in the case of GRB- 
DLA 050820 but predominantly negative velocities for 
GRB-DLA 051111. Therefore, one example exhibits mo- 
tions reflective of a random (i.e. virialized) velocity field 
and the other inspires an outflow interpretation (see also 
Thoene et al. 2006). It is too early to draw generic 
conclusions about the velocity fields of GRB-DLAs, but 
future observations akin to Figure 9 will offer an unprece- 
dented view into the nature of gas kinematics at high z. 

Let us now address the offset in the W1526 distributions 
between the GRB-DLAs and QSO-DLAs (Figure 4b). 
Could the offset be explained by differences in sightlinc 
configurations? As with the Av 9 o values, there are com- 
peting effects. The bias to smaller impact parameter will 
only be important if the Si + velocity fields have signifi- 
cant variations on kpc scales. This could occur if the ma- 
jority of gas is located near star-forming regions, i.e. the 
origin of the GRB sightline. If the Si + gas is distributed 
throughout the halo (i.e. tens of kpc), then an impact 
parameter bias should be unimportant. Meanwhile, if 
the velocity field has significant contribution from in- 
fall or outflow then restricting the GRB-DLA sightlincs 
to originate from within the ISM leads to a systematic 
underestimate. Therefore, sightline geometry will bias 
the W1526 statistic to larger values for the GRB-DLA 
only if the velocity field has significant variations on kpc 
scales and is asymmetric. Numerical simulations do in- 
dicate gas being accreted tends to fall in along radial 
orbits (A. Dckcl, priv. comm.). This may imply larger 
11^1526 values for sightlines originating at the center of the 
galaxy (GRB-DLAs) as opposed to large impact param- 
eter (QSO-DLAs). While differences in sightline config- 
urations may contribute to the differences in W1526 be- 
tween QSO-DLAs and GRB-DLAs, we contend that the 



observed differences in the W\$2Q distributions for GRB- 
DLAs and QSO-DLAs reflect actual differences in the 
gas kinematics of their host galaxies. 

A possible explanation for larger VF1526 values in GRB- 
DLA is that their host galaxies are significantly more 
massive. There are two arguments against this asser- 
tion. First, the Av 9 o distributions of the QSO-DLAs 
and GRB-DLAs are similar (Figure 4a) . While the T / Fi526 
statistic may be a better tracer of the gravitational po- 
tential (see below), the Av 90 statistic should also dis- 
tinguish between large differences in dark matter halo 
masses. Second, direct imaging of GRB host galaxies 
reveal they have low luminosity (Le Floc'h et al. 2003; 
Fruchtcr et al. 2006) and, presumably, low mass. We 
consider it unlikely that the QSO-DLA have substan- 
tially lower luminosity (and mass) than that observed 
for the GRB-DLAs. 

The most likely interpretation for the higher 14^1526 
values in GRB-DLAs is that the kinematics of gas con- 
tributing to the W1526 statistic has significant contribu- 
tion from galactic-scale outflows. We are guided to this 
conclusion in part because the GRB phenomenon is di- 
rectly linked to active star formation, i.e., GRB host 
galaxies have systematically higher current star forma- 
tion rates than QSO-DLA galaxies. One observes that 
GRB host galaxies have star formation rates which are 
large for their luminosity (Christenscn, Hjorth & Goros- 
abel 2004) . This characteristic is described as a high spe- 
cific SFR and it separates GRB host galaxies from the 
general population of high z galaxies. Therefore, a sce- 
nario where processes related to star formation stir the 
velocity fields traced by Si + gas may naturally account 
for the systematic offset in between QSO-DLAs 

and GRB-DLAs. These processes may include galactic 
outflows inspired by major mergers (Cox et al. 2006) or 
galactic fountains driven by supernovae feedback (Mac 
Low & Ferrara 1999; Dong, Lin & Murray 2003). The 
latter effect may be challenged, however, by the short 
lifetime expected for GRB progenitors. In this respect, 
the velocity fields of GRB-DLAs may more resemble the 
outflows of the Lyman break galaxies (e.g. Pettini et al. 
2002). We eagerly await new observations like those pre- 
sented in Figure 9 to further test these scenarios. 

5.3. Interpreting the Kinematic- Metallicity Correlations 

in DLAs 

In Section 4.2 (Figure 6), we presented the observed 
correlations between the kinematic statistics (Av 9 o, 
W 1526 ) and the ISM metallicity [M/H]. The Av 90 -[M/H] 
correlation has been previously identified for QSO-DLA 
samples (Wolfe & Prochaska 1998; Ledoux et al. 2006). 
Although its scatter is substantial, the statistical signif- 
icance is high and Ledoux et al. (2006) report a best- 
fit power-law of the form [M/H]cx 1.51ogAw 90 . These 
authors have interpreted the result in terms of a mass- 
mctallicity relation, i.e. the Av 90 statistic traces the grav- 
itational potential of the DLA galaxy. Our results on the 
GRB-DLA lend further support to this interpretation. In 
particular, we have shown that the At)go statistic is dom- 
inated by material near the GRB, i.e. the neutral ISM. 
Although ISM velocity fields are stirred by supernovae, 
they are generally dominated by gravitational motions 
(e.g. rotation, turbulence, accretion). 

A key aspect of the Au 90 -[M/H] trend is its large scat- 
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ter which cannot be attributed to observational uncer- 
tainty. We interpret the scatter as a natural consequence 
of the QAL experiment. If the velocity field is dominated 
by rotation, then a large scatter results simply from a 
range of impact parameters and disk inclinations (e.g. 
Prochaska & Wolfe 1997). In this case, only the upper 
bound of the Avgo,[M/H] locus would describe the av- 
erage rotation speed of galaxies at a given metallicity. 
Similarly, if the velocity field results from turbulence or 
accretion, then At>go will be sensitive to the number of 
clouds intersected by the sightline. In this scenario, how- 
ever, one would require large variations in the number of 
'clouds' intersected by sightlincs, which is not entirely 
supported by the observations (e.g. Prochaska & Wolfe 
1997). 

In contrast to the Au 90 -[M/H] trend, the W 1526 ,[M/H] 
pairs for the QSO-DLAs follow an extremely tight cor- 
relation (Figure 6b). In fact, most of the scatter may 
be attributed to observational uncertainty. In several 
respects, this is a stunning result. The [M/H] values 
represent the average metallicity along the entire sight- 
line, which is dominated, of course, by the gas corre- 
sponding to the majority of the total optical depth, i.e., 
the gas which defines the Aw 90 statistic. In contrast, 
gas which has negligible contribution to the measured 
metallicity can dominate the W1526 statistic (Figure 9). 
Furthermore, we have shown (for the GRB-DLAs) that 
[M /H] is set by gas local to the ISM of the galaxy whereas 
the IF1526 statistic has significant contribution from gas 
external to the ISM. Indeed, quasar absorption lines 
are frequently associated with halo gas surrounding a 
galaxy (e.g. Bergeron & Boisse 1991; Steidel 1993; Chen, 
Lanzetta & Webb 2001; Charlton & Churchill 1998). 
Therefore, in the DLAs we reach a rather surprising con- 
clusion that the local ISM properties (metallicity) are 
tightly correlated with the large-scale velocity field. 

Let us now comment on the origin of the Wi526-[M/H] 
trend. The first point to stress is that we derive a simi- 
lar power-law as observed for the Ai>go,[M/H] pairs in 
QSO-DLAs (see also Murphy et al. 2007). We con- 
clude that the kinematic-metallicity trends have the same 
physical origin. One may gain insight by comparing 
our trend with the velocity-metallicity correlations ob- 
served for other galactic populations. Dekel & Woo 
(2003) have found the velocity dispersion V in dwarf and 
low-luminosity low surface-brightness galaxies at z = 
correlates with the stellar mass M„: V oc Mf 24 . Fur- 
thermore, their compilation shows a tight correlation be- 
tween metallicity and stellar mass, M/H oc M, im- 
plying a metallicity-velocity correlation: M/H cx V' 1 " 66 . 
The correspondence between the kinematic-metallicity 
trend of local low-mass galaxies with that of the QSO- 
DLAs is remarkable. Although one refers to this trend 
locally as a mass-metallicity relation, it is important to 
stress that simple scaling laws do not predict the ob- 
served trend. Dekel & Woo (2003) have shown that 
in a closed-box model the metallicity is independent of 
stellar mass; indeed, this follows observations for mas- 
sive, high surface brightness galaxies (e.g. Tremonti et 
al. 2004). Therefore, Dekel & Woo (2003) argued that 
the mass-metallicity trend observed for lower mass galax- 
ies requires supernovae feedback (but see also Lee et al. 
2006; Tassis, Kravtsov & Gnedin 2006). In this respect, 
therefore, one may speculate whether the velocity fields 



observed in the QSO-DLAs represent gravitational dy- 
namics or feedback processes. Let us now consider argu- 
ments for and against each interpretation. 

There are several arguments supporting the interpreta- 
tion of Figure 6 in terms of a mass-metallicity correlation, 
i.e. that the velocity fields trace the gravitational poten- 
tial. First, the observed slopes of the Aiigo and W1526- 
[M/H] trends follow the same trends as local dwarf and 
low-luminosity low-surface brightness galaxies (Dekel & 
Woo 2003) whose velocity dispersions presumably reflect 
their gravitational potentials. Second, the large scat- 
ter in the Awgo-[M/H] trend can be explained by varia- 
tions in impact parameter and inclination through the 
DLA galaxies. Third, a tight trend is reasonable for 
Wi526>[M/H] pairs if the VF1526 statistic is dominated by 
halo dynamics. For virialized motions characterized by 

1 /2 

a velocity dispersion a viria i, Wi 52 6 oc nj ld a viria i and the 
scatter in this quantity will be small for a given velocity 
dispersion <J V i r iai if the number of clouds penetrated by 
the sightline n c id 3> 1. The latter point is supported by 
the observations (e.g. Figure 1). Fourth, the Si II 1526 
profiles of at least some GRB-DLAs show both negative 
and positive velocities relative to the systemic (Figure 9) . 
This is expected for gas that traces a virialized (random) 
velocity field (see also Robertson & Chen, in prep.). Fi- 
nally, Moller, Fynbo & Fall (2004) have noted that the 
only DLAs that have been successfully imaged at z > 2 
also have high metallicity, suggesting stellar mass may 
correlate with enrichment. 

There are two counterarguments to interpreting the 
Av 9 o and W1526 statistics as tracers of the gravitational 
potential of individual dark matter halos. First, it is may 
be unrealistic to explain velocity widths that exceed sev- 
eral hundred km s _1 in terms of a single galactic poten- 
tial. At very large Augo or W1526 values, therefore, we 
expect that additional velocity fields contribute. These 
could include galactic-scale outflows but also the pecu- 
liar motions of multiple galaxies along the sightline (e.g. 
Mailer et al. 2001). Second, Bouche et al. (2007) have 
reported an anti-correlation between Mg II equivalent 
width 9 and the cluster-length of these absorbers (a subset 
include DLAs Prochter et al. 2006) with large, red galax- 
ies (LRGs) . These authors have interpreted this result in 
terms of a wind scenario: systems exhibiting larger equiv- 
alent widths occur in less massive, star-bursting galaxies. 
We note, however, several caveats to their interpretation: 
(i) the Mg II absorbers have z < 1 and may have very 
different characteristics (star-formation rates, metallic- 
ity, gas content) than the z > 2 QSO-DLAs that are 
plotted in Figure 6; (ii) Zibetti et al. (2007) find the 
integrated light of galaxies corresponding to Mg II ab- 
sorbers is brighter for larger equivalent width systems 
suggesting larger masses. 

Now consider the evidence in favor of interpreting the 
VF1526 statistic (and perhaps Av 90 ) as galactic-scale out- 
flows. First, as noted above, the trend of decreasing 
correlation-length with increasing equivalent width offers 
at least qualitative evidence for a wind scenario (Bouche 

9 In absorbers where one observes both the Si II 1526 transition 
and the Mg II doublet, one finds the line profiles track each other 
very closely. As such, their equivalent widths are highly correlated 
and roughly scale as the inverse of their wavelengths: W1526 ~ 
VK2796A1526/A2796 (in detail we find W2796 ~ 3Wi 5 26)- 
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et al. 2007). Second, the kinematics of some strong 
Mg II absorbers are consistent with that expected for 
symmetric outflows (Bond et al. 2001). Third, if star 
formation drives galactic-scale outflows, then one may 
expect a correlation between metallicity and Wi526- It 
may be difficult to precisely derive the observed trend 
([M/H] oc 1.5 log W1526) from first principles, but we can- 
not rule out such a scenario. Third, we find the Lyman 
break galaxy MS 1512-cB58 (Pettini et al. 2002) lies di- 
rectly on the Wi526-[M/H] trend of the QSO-DLAs (black 
diamond in Figure 6b). Given that the LBG gas kine- 
matics are dominated by outflows (Pettini et al. 2002), 
the correspondence between the LBG data and the QSO- 
DLA trend is at least suggestive that outflows may play 
a role in shaping the observed correlation. The com- 
parison with the QSO-DLAs may not be appropriate, 
however, because the LBG sightline is restricted to orig- 
inate at the central, star-forming region of the galaxy. 
This implies a special impact parameter and one that 
presumably probes only half of the total velocity field. 
In these respects, a comparison with the GRB-DLAs is 
more appropriate and the correspondence with the QSO- 
DLA trend may be simple coincidence. 

There are several arguments against interpreting the 
W1526 measurements as galactic-scale outflows. First, 
the similarity in the W1526 and Ai>go-[M/H] trends would 
require one to interpret both statistics in terms of out- 
flows (or invoke coincidence). In turn, one would have to 
argue that these winds have gas with large neutral frac- 
tions (Prochaska et al. 2002) and significant cross-section 
to gas with N m > 10 21 cm~ 2 (e.g. Schaye 2001). Sec- 
ond, we believe a galactic-scale outflow scenario would 
be severely challenged to reproduce the tight correlation 
observed for W1526 and metallicity. Assuming outflows 
are generally collimated, the VF1526 statistic should ex- 
hibit a wide range of values for a given wind speed as 
a function of the sightline configuration (i.e. impact pa- 
rameter and inclination). Furthermore, the W1526 value 
is given by the internal dynamics of the wind (i.e. its tur- 
bulence) because it is a relative measure independent of 
the systemic velocity. It is not obvious that this quantity 
could be tightly correlated with metallicity. Third, there 
are at least some GRB-DLAs where the W1526 velocity 
field is not strictly consistent with outflows (Figure 9). 
Fourth, mass-metallicity relations have been identified 



for numerous galaxy populations at a range of redshifts 
(Dekel & Woo 2003; Tremonti et al. 2004; Kobulnicky & 
Kewley 2004; Savaglio et al. 2005; Erb et al. 2006). It 
would be very remarkable for the QSO-DLAs to exhibit 
a tight velocity-metallicity trend that was unrelated to 
a mass-metallicity correlation. Fifth, we question a sce- 
nario which requires galactic-scale outflows in all galaxies 
tuned to give both low and large W1526 values with such 
small scatter. Sixth, comparisons of Mg II absorption 
line profiles with H II emission-line kinematics are in- 
consistent with a wind scenario and support dynamical 
motions (Steidel et al. 2002). Finally, while observers 
are quick to invoke outflows to explain extreme veloc- 
ity fields, we question whether outflows can reasonably 
explain widths approaching and in excess of 1000 km s" 1 . 

On balance, we currently favor the interpretation of 
the velocity fields in QSO-DLAs as being dominated by 
gravitational motions, at least for W1526 < 1.5A and 
Avgo < 200 km s _1 . At large VF1526 and large Awg , it is 
reasonable to expect that additional velocity fields con- 
tribute. Indeed, an inspection of Figure 6b reveals that 
the QSO-DLA galaxies with largest departure from the 
[M/H]-Wi526 relation occur at large W1526 values. These 
few QSO-DLAs have especially large VF1526 for their ob- 
served metallicity and/or low metallicity for the observed 
W1526 value. We suspect that the former inference is 
the correct interpretation and that these systems have 
contributions to the W1526 statistic beyond the halo gas 
dynamics. These QSO-DLAs are the most viable can- 
didates for sightlines penetrating galactic-scale outflows. 
This interpretation could also explain the frequency of 
GRB-DLAs that lie off the QSO-DLA trend. As noted 
in the previous section, GRB-DLAs have especially large 
specific star formation rates. Therefore, GRB may flag 
galaxies at a time when they are most likely to drive 
galactic-scale outflows. Given the identification of GRB 
with massive, short-lived stars, there may not be suffi- 
cient time to drive these winds with supernovae feedback. 
Instead, one may need to consider merger-driven, galac- 
tic outflows. 
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